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Abstract. We study the dynamics of a quantum particle in R"+'" constrained by 
^ . a strong potential force to stay within a distance of order h (in suitable units) from 

I a smooth n— dimensional submanifold M. We prove that in the semiclassical limit 

' the evolution of the wave function is approximated in norm, up to terms of order 

, h}/"^, by the evolution of a semiclassical wave packet centred on the trajectory of the 

^■f-^ ' corresponding classical constrained system. 

o ■ 

r-| ■ 1- Introduction 

> 

<^ 

The aim of this paper is to study the semiclassical limit of a nonrelativistic quantum 
Hamiltonian system in the configuration space M"+"^, constrained to a submanifold 
^ ■ M C by a confining potential which becomes infinite, in a suitable sense to be 

' defined, when we move away from M. 



We derive an effective Hamiltonian for the classical motion on M, using the 
technique developed in a series of papers by Hagedorn (Hagedorn 1994, Hagedorn 1998 
and references therein) to construct approximate solutions to the Schrodinger equation 
which are localized along a classical trajectory. 

We limit ourselves to Hamiltonians of the form 

if, = M! + Viq) + W%q) 

p := —ihVq q := q- 

where is the confining potential and e is a small parameter which we will make 
eventually go to zero (in section 01 we examine the motion of a particle in a magnetic 
field, which, under suitable conditions, can be put in form (jl.lj) ). 

To explain the characteristic features of the method we employ, we first analyze in 
detail a number of explicit cases (M" embedded into R*^"*""*, a smooth curve embedded 
into a plane); we show then how the procedure generalizes to (non-fiat) submanifolds of 
arbitrary dimension and codimension. 
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The confining potential approach to imposing a constraint has been used often in 
the hterature for a variety of reasons. 

In classical mechanics, it has been employed mostly to "realize holonomic 
constraints" (Frocsc and Herbst 2001), i . e., to justify the use of D'Alembert principle 
in deriving the Lagrange function for systems subject to holonomic time-independent 
constraints (which was the starting point for the research performed in Takens (1980)). 
Other traditional applications include the analysis of magnetic traps and mirrors, whose 
first complete mathematical discussion was given in Rubin and Ungar (1957), which was 
also the first rigorous investigation in the field (a detailed treatment of these problems 
from the point of view of weak convergence, with extensions to arbitrary Riemannian 
manifolds and molecular dynamics, can be found in Bornemann (1998)). 

In quantum mechanics, the limit of large restoring force has been considered 
essentially for two reasons. 

The first is that it offers a way, different from the intrinsic one (Henneaux and 
Teitelboim (1992) and references therein), to quantize constrained systems (da Costa 
1981, da Costa 1982, Jensen and Koppe 1971, Kaplan et al 1997, Maraner and Destri 
1993, Mitchell 2001, Schuster and Jaffe 2003). 

The second is that in mesoscopic physics {i. e., the branch of physics which studies 
small objects, like thin films and quantum wires) there exist systems which have one, or 
more, dimensions much smaller than the others, and arc then well described, in a zero 
order approximation, by an {n — A';)-dimcnsional confined system, 1 < k < n — 1, (for 
the physical background and mathematical models see Duclos and Exner (1995), Exner 
(2003) and references therein). 

It was only recently (Proese and Herbst 2001, Teufel 2003) that a comparison 
between the classical and the quantum case was attempted. 

The main problem one runs into is that, due to the Heisenberg principle, the mean 
value of the Hamiltonian operator diverges for every initial condition in the constraining 
limit (the better we localize the wave function on the submanifold M, the bigger the 
mean of the square of the momentum becomes), while almost all theorems available in 
classical mechanics deal with finite energies. 

To overcome these difficulties, Frocsc and Herbst state and prove a theorem on the 
classical case with unbounded energy, which, however, does not seem very natural from 
a physical point of view, while in Teufel (2003) it is suggested to consider, instead of the 
limit of large restoring forces, the limit of weak forces in the non-constraining directions 
(they are equivalent in classical mechanics, up to a rescaling of space-time). 

We propose a different approach, based on the fact that in quantum mechanics 
there exists an a priori length scale defined through h (in units in which time and mass 
are of order one). 

In real systems, like the mesoscopic ones mentioned above, the transversal directions 
contain at least some atoms, so any realistic layer cannot become smaller than h, which 
is the order of magnitude of atomic dimensions. Therefore, in our opinion, it is necessary 
to link the squeezing scale, determined by the constraining potential, to the quantum 
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scale given by h. 



1.1. A simple example 

To illustrate this point, we consider the standard two-dimensional example 

Hs = + ^-(x) V, (1.2) 

where : R ^ is an arbitrary smooth function which satisfies uj{x) > u;* > 

Vx e M. 

The squeezing scale is determined by e, and we want it to be a function of h, 
e = e{h). Since, as we argued before, e cannot become smaller than h, and it has to go 
to zero when h ^ Q (to achieve the constraining limit), the simplest choice is 

e = aK' 0<a<l a fixed > (1.3) 

(there is no loss of generality, since what matters is the behaviour of e{fi) when ^ ^ 0). 
With this choice, the Hamiltonian p.2|) becomes 

and we want to examine the limiting behaviour of the dynamics generated by Hfi, when 
If we unitarily scale the transversal direction to factorize h 

we get 
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5» - -\dl + fi- 



If a 7^ 1, using the same techniques illustrated in the next sections, it can be shown 
that the infiuence of the normal motion on the longitudinal one is suppressed, and the 
effective Hamiltonian is the free one. Therefore, in the following, we consider only the 
more interesting case a = 1. 



1.2. Outline of the paper and summary of results 

In the next section we analyze a generalization of (jl.2j) . studying the case of a potential 
confining to a fiat sub manifold M of M""*"*". We realize the constraining limit through 
dilations in the direction normal to M, i. e., we put W^{x,y) = W{x,y/e). This allows 
us to consider generic dependence on the transversal variables, unlike what is usually 
made in the literature (Bornemann 1998, Froese and Herbst 2001, Takens 1980), where 
the first non zero term in the Taylor expansion of the potential around the constraint is 
the quadratic one, and so the problem is reduced to the analysis of harmonic motions. 

In section El we consider a two-dimensional example where the constraining limit 
is realized through the more traditional method of scaling of the coupling constant, i. 
e., W'^{x,y) = e~'^W{x,y). In the case of a spectrally smooth potential confining to 
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a nondegenerate critical curve (for the definitions, see Appendix A and Appendix B ) 



the semiclassical limit motion we get along M is the same as the homogenized classical 
motion found by Bornemann (1998). 

In section |31 we show that an analogous result holds for an n-dimensional 
nondegenerate critical submanifold embedded into R*^"*"™. We exploit Hagedorn's 
multiple scale technique to construct squeezed states whose centre and dispersion take 
account of the (non-trivial) curved background. 

Finally, we address an Hamiltonian showing the phenomenon of Takens chaos 
(Bornemann 1998, Takens 1980), which is encountered when the constraining potential 
is not spectrally smooth. In classical mechanics, the motion on the submanifold M is 
not deterministic anymore, i. e., it is not described by a natural mechanical system 
on M and the limit set obtained forms a funnel. We show that the semiclassical limit 
offers a natural way to reduce (but however not to eliminate, in general) the degeneracy, 
linking different trajectories in the funnel to different quantum initial conditions. 

2. Constraints by normal dilations 

Let M = W and W%q) = W{x,y/e), where we split q E as {x,y), x G M", 

y e M"". 

We suppose that 

V,W E Lf^^ and are bounded from below, (2.1a) 
lim W{x,y) = oo Vx G M"* (confining hypothesis). (2.1^) 

We impose also an implicit smoothness hypothesis on the potentials, through a 
condition on the resolvent of the reduced Hamiltonian i){x), to be defined below. 

As argued above, we put e = ah. Actually, since we have several normal directions, 
we can choose different e/h ratios for each one. 

Defining 

i/-=f!r ••• ^) (2-2) 



equation ()3.H1 becomes 

^ft = ^ + V{x) + W{x, ya/h), (2.3) 

where, for the sake of simplicity, we suppose that V{q) does not depend on y. 
Scaling the transversal directions by the dilation operator 

(D,V^)(x,y)=7'"/V(a:,71/), (2.4) 

we get an Hamiltonian of the same form as the Born-Oppenheimer operator, used in 
molecular physics, 

D[_,HnDn-^ =: Hbo = -^^x + f)(x), 

1 2 (2.5) 

[){x) = --Ay + W{x,ya) + V{x). 
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It follows from ipTaj) . (|2?n|) that i){x) is, for each well defined self-adjoint 
operator, with compact resolvent and nondegenerate ground state (Reed and Simon 
1978). 

We suppose in addition that has a smooth dependence on x, namely that 
— i)~^ is a function of x, for some / > 2. This makes its eigenvalues E(x) 
(which we will call also "transversal" or "normal" energy levels) C' functions of x away 
from crossings or absorption in the continuous spectrum. 

The behaviour of Born-Oppenheimer Hamiltonian when h is well understood 
(Hagedorn 1994, Teufel 2003). 

The transversal motion adiabatically decouples from the longitudinal one and stays 
approximately in a bound state of i){x) for a fixed value of x. On the other hand, the 
longitudinal motion depends on the transversal one because it feels an effective potential 
which is equal to the normal energy. 

Using standard results (Hagedorn 1994) we can elaborate on this qualitative 
argument: 

Theorem 2.1. Suppose that there exists an open set f/ C such that l){x) has 
a nondegenerate eigenvalue E{x) for x E U, with corresponding real normalized 
eigenfunction $(a;). 

Let a{t) and r]{t) be the solutions of the classical equations of motion with potential 
E{x) (which exist and are unique since E{x) is C\U) and hounded from below) 

a{t) = r^it) (2.6) 

m = -VE{a{t)), (2.7) 

a(0) = ao v{0) = Vo, (2-8) 

then, for t G [0,T], 

exp ( - '-^Hn^MMO), B{0), h, a(0), 7^(0), x)F(x)Da-i$(x) 

- exp U^^^^^iAit), B{t), h, a{t),r]{t),x)F{x)Dn-Mx) 

= 0{h}/^), (2.9) 

where S{t) is the classical action, A{t) and B{t) are linked to the dispersions of (fk 
(respectively) position and momentum and F is a cut function which is zero outside a 
neighbourhood of the classical trajectory {a{t) : t G [0,T]}. 

Remark 2.1. The functions (pk{^,B,h,a,ri,x) were introduced by Hagedorn, to whom 
we refer for the notation (Hagedorn 1998). They are a useful tool in studying the 
semiclassical limit of quantum mechanics and they coincide with the "squeezed states" 
widely used in quantum optics (Combescure 1992). Essentially, they are minimal 
uncertainty wave packets with different spreads in position and momentum. 

Remark 2.2. We will give a proof of a slightly more general version of theorem 1)2.11) in 
sections El and m where we analyze the Laplace-Beltrami operator in a curved space. 
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2.1. Comments and examples 

Let us analyze in greater detail the approximate evolution found in ()2.9|1 . 

The transversal wave function Dfi-i^{x) clearly describes a motion confined to the 
submanifold M = R" , since 

< y >=< Dn-i<^{x),yDn-iHx) >= h < <l>(x),?/$(x) >= O(^) 

(A^,)2 =< Dn-iHx),y^Dn-i^x) >-< Dn^i^x),yiD^-.^x) >^= 0{h''), 

while both < py > and < Apy > are 0(1). 

One should note however that we did not require W to have a strict minimum on 
M. Actually this is not needed, since in our scale the average position of the normal 
motion is always "seen" to be approximately zero, as equation fl2.10|) shows. 

In the standard case where 

W{x,y) = -J2^,{xryl (2.11) 

j=i 

the effective potential for the motion on M will be 



(2.12) 



This is exactly the homogenized potential found by Bornemann (1998) and Takens 
(1980), where the are, in the classical case, the adiabatic invariants associated to the 
normal oscillations {i. e., the energy-frequency ratios). 

Varying the squeezing factors a^, or the transversal wave function $(x), i)i can be 
made to assume every positive value (the value 'di = can be obtained suppressing 
the ith mode as we explained in section [T!T|) . The harmonic potential is particular in 
this respect, because, as far as the effective potential is concerned, all normal states are 
equivalent, since the various choices for $(a;) correspond simply to suitable scalings of 
e and h. 

One could even use an x-dependent scale, e = a{x)h, without altering substantially 
the structure of equation ()2.12|) . 

Such a simple picture cannot be expected when W is not harmonic. 

In general, the effective potential will have a non-trivial dependence both on the 
parameters a := (ai, . . . ,am) and the transversal wave function. This gives a host of 
well-defined classical motions on M, whose form, however, cannot be given explicitly as 
in the harmonic case. 

It would be interesting, for instance, to compare the semiclassical effective 
Hamiltonians produced by a "flat" confining potential, like the sextic harmonic 
oscillator, 

W{x,y) = V^{x)y^ + VQ{x)y^ {x,y) eR^ Ve{x)>V.>0 (2.13) 
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with the corresponding homogenized classical motions (if any exists), to see if it is 
possible to reproduce them in a purely classical way. 

Unfortunately, the spectrum of the reduced Hamiltonian associated to ()2.1H|1 is 
known only for particular values of the squeezing parameter a. For example, if a = 1 
and V4(a;)^ = 12F6(x)^/^ it is known (Skala et al 1996, Ushveridze 1994) that the ground 



but it is not possible to write an explicit expression for all values of a. 

3. Constraints by scaling of coupling constant: a curve in a plane 

In this section we analyze, in a fairly detailed way, a two-dimensional example where 
= e~^W . It allows to explain the main differences between the curved and the flat 

case, avoiding technical complications arising from higher codimensions, which are not 

essential for the result, and will be illustrated in next section. 

We suppose, in the same spirit of ()2.1a|) . that V and W are C°° and non-negative, 

but, as is customary in classical mechanics (Bornemann 1998, Takens 1980), we replace 

fl2.16|) with the hypothesis that is a spectrally smooth potential constraining to a 

nondegenerate critical curve M( |Appendix A| and [Appendix B| ). 

Our starting Hamiltonian (with the prescription e = ah) will be then 



Squeezed states are particularly suited to studying this sort of situations, where M 
is not fiat, because, as ()2.9p shows, the evolution of a localized state is approximately 
described (for a bounded time interval) by localized states. This allows us to analyze the 
motion using one coordinate chart only and therefore local expressions for the operators 
involved. 

Essentially, what we will do here is to adapt the arguments of the last section to a 
curved case, constructing an approximate solution to the Schrodinger equation which, 
in suitable coordinates, is still given by a squezeed state in the longitudinal direction 
and an (harmonic) oscillation in the transversal one. 

3.1. The Hamiltonian in curvilinear coordinates 

We fix a tubular neighbourhood of M, and we consider a single chart of tubular 
coordinates, defined on C 1^. 

This simply means that, given a local parametric representation of M in terms of 
its arc length s, qui^s) = (xa/ (s), ?/a/(s)), we can write (for g G '^) 



state is 




(2.14) 




(3.1) 




(3.2) 



where ti{s) is the unit normal of M. 
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In writing (|3.2p we used the natural linear structure of tubular coordinates. A more 
invariant, but less manageable, relation would be 

q = expg^^ qM G M G Tg^,M^, 

where exp is the geodesic exponential map (Lang 1995). In the following, however, we 
will stick to (jSI2I). 

When q varies over s and u vary, respectively, over two intervals / and J. 
Lemma 3.1. The Hilbert space L'^{'^,dq) is isometric to L^{I x J,dsdu). 

Proof. This well-known lemma results from two facts. 

First, the choice of curvilinear coordinates provides an isometry of L^{'^,dq) to 
L^(/ X J, g^/'^dsdu), where 

(7^/2 = 1 - k{s)u (3.3) 

is the Jacobian of the transformation {x,y) — > {s,u), and k{s) is the curvature of M. 

Second, the multiplication by g^^^ is a unitary operator from L^(/ x J, g^^'^dsdu) to 
L^(/ X J, dsdu). □ 

In the following, we will denote the isometry constructed above by U : L^(^, dq) — > 
L^(/ X J, dsdu). 

We remark that U maps C^i"^) onto C^{IxJ) and Hn maps C^i"^) into C^{'^), 
so, denoting, with abuse of notation, the restrictions of U and Hn to functions with 
the same symbols, we have 

UHnU^ : C^{I X J) ^ C^{I X J) 



2 (1-A;(s)m)V2 ^Vl-fc(s)M '(1-A;(s)m)i/2 

- ^TT^^TTV^ + Vis, u) + {ah)-^W{s, u) (3.4) 



2 " 8 {l-kis)u] 



r5.^-77^^^^5.-^^^?(^,-) (3.5) 



2{l -k{s)uy ' {l-k{s)uf 

'Y^l + V{s, u) + {ahr^W{s, u), (3.6) 

where V and W are and written in curvilinear coordinates and h^Q is an 
extrapotential of purely quantum origin which depends on the curvature k{s) (da Costa 
1981, Jensen and Koppe 1971). It appears also in mesoscopic physics, and can give 
rise to interesting phenomena, hke bound states, in a quantum waveguide (Duclos and 
Exner 1995). However, it will not concern us, since it disappears in the lowest order of 
semiclassical approximation. 

Using again a dilation operator in the transversal direction u, 

: L'^{I X J^, dsdu) L?{I x J, dsdu) 
{D^iIj){s,u) = 7^/^(5, 7m) 

Jy := {'ju : u G J}, (3.7) 
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we get the final form of the Hamiltonian which we will employ in the estimates: 
Hbo : C^{I x J^-i) ^ C^{I x J^-i) 



Hbo -.^Dl.UHnU^D^-, 



u 



9; - TI-^TtV^^. - ^ Qi^^ H + (3.8) 



2(1 -hk(s)uy ' (l-hk(s)u)^ 
where 

[){s) = ~dl + V{s, hu) + {aIi)-^W{s, hu). (3.9) 
Remark 3.1. Note that 

1 ~ h ~ T rhu 

iah)-'Wis, hu) = ^dlWis, 0)u' + ^^X^, 0)u' + ^ dv{hu - vfdiW{s, v) 

= ^^(«) V + -^dlW{s, 0)^3 + R,^h, u), (3.10) 

V{s, hu) = V{s, 0) + hudj'is, 0) + / dv{hu - v)dlV{s, v) 

Jo 

= V{s, 0) + n^/9„y(s, 0) + Ri{h, u). (3.11) 

The scaling in the normal direction eliminates the dependence of f) on only at the 
lowest order in the Taylor expansion around the constraint (which is the quadratic one 
since M is a nondegenerate critical curve). 

Prow now on, we will denote by i)^'^\s) the harmonic part of i){s): 

F^)(.) := ~dl + ^o;(.) V + V{s, 0). (3.12) 

3.2. The approximate evolution 
In this subsection wc prove the 

Theorem 3.1. Let $(s, -u) be a real normalized eigenstate of l)^'^\s), considered as 
an operator on L'^{W,du), with eigenvalue E{s). Let a{t) and rjit) be the solutions of 
the classical equations of m,otion with potential E{s), and let F{s,v) be a function in 
C^{I X J) which is equal to 1 for s in a neighbourhood of the trajectory {a{t) : t E [0, T]} 
and V near to 0. 
Then 

exp ( - '-^Hn)U^Dn-iMMO), B{0), h, a(0), r/(0), s)F{s, hu)<^{s, u) 



h 



- exp U^Dn-i^k(A(t), B(t), h, a(t),ri{t), s)F(s, hu)^s, u 

= 0(n^/2), (3.13) 
where S{t) is the classical action associated to {a(t),ri{t)) . 

Remark 3.2. The function ipkiA{t), B{t), h, a{t),r]{t), s)F{s, hu)^{s,u) is in C^{I x 
Jft-i), so UWfi-icpk{A{t),B{t),h,a{t),7]{t),s)F{s,hu)^{s,u) belongs to C^{'^). 
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The proof will follow closely the pattern developed by Hagedorn (Hagedorn 1994), 
but the remainder we get is different from that found by him, since f) contains terms of 
order h and the kinetic part of ()3.8|1 is not simply —(Ji?/2)dl. 

The basic tool we use is a simple application of the fundamental theorem of calculus 
(also known as Duhamel formula). We give it without proof. 

Lemma 3.2. Suppose is a family of self-adjoint operators for h> 0. Suppose ip{h, t) 
belongs to the domain of Hf^, is continuously differentiate in t, and approximately solves 
the Schrddinger equation in the sense that 

ihdtip{h, t) = Hni^ih, t) + C{h, t), (3.14) 

where ({h,t) satisfies 

\\C{hM<K^t) (3.15) 

for <t <T. Suppose \E'(^, t) is the exact solution to the equation 

ihd,^{h,t) = Hn^{h,t) (3.16) 

with initial condition \E'(^, 0) = ipih.Q). 
Then, for < t < T , we have 



|^(n,t) -7/'(n,t)|| < r 



drfi{h, t). 



(3.17) 



Suppose now that iljap{s,u,t) G C^{I x J^^-i) is an approximate solution to the 
Schrodinger equation associated to the local Hamiltonian ()3.8|) . 



with 



for 



< t < T. 



(3.18) 
(3.19) 



L2 



1/2^ 



This implies that 

ihdfi^Dn-^i,^p = HnU^Dn-i^p,^ + ah,t), 

with||c(n,t)|U2(^,) = o(;i3/2)_ 

Using lemma ()3.2|) we get finally 

I exp ( - '^H,^U^Dn-i^^p{t = 0) - f/^Da-i^ap(t) 

Therefore, to prove theorem ()3.1|) we will construct an approximate solution to 
dSIHl), of the form 

V'ap(s, M, t) = iJo{s, u, t) + h^j^is, u, t), (3.20) 

with ipo{s, u, t) = exp{iS (t) / h)ipk{A{t) , B(t), h, a(t),ri(t), s)F{s, hu)^{s, u) (the notation 
t/'^ means that the transversal part of this term is orthogonal to $). 

An educated guess about the form of the remainder t/)^ can be made employing a 
multiple scale technique, which allows to split the adiabatic and the semiclassical scale. 



dlW{s,0)u^^{s,u)], (3.21) 
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We will elaborate on this procedure in the more complicated case of next section, 
so here we limit ourselves to verify that the right choice is 

^2^(s, u, t) = ipk{A{t), B{t), h, a{t),r]{t), s)F{s, hu) 

xf{s) 

1 

where f(s) is the bounded inverse of the restriction of [[^'•^''(s) — E(s)] to the orthogonal 
complement of in L'^{W,du). 

Estimate ()3.19j) will follow if we note the following facts: 

(i) the terms containing derivatives of F are 0{h°°). For instance, 

I dsdu\d^F{s,hu)<^k{s)d^'l>{s,u)\^ (3.22) 

= I dsdv\d^F{s,v)yDk{s)h^/^d^<^>{s,vh-^)\^ <exp{-Ch'^), 

JlxJ 

since d^^F has support away from zero in v, and 9^$(s, is a polinomial times 

a Gaussian, in u = vh~^. 

The derivatives with respect to s can be estimated in the same way, since ipk is a 
Gaussian in [s — a(t)]//i^/^. 

(ii) The term 

(T^^M^ ' 

is 0(^^) since d^^pk is 0{h^^). 

(iii) The term 

is 0(^^) since Q{s,u) is bounded on the support of F. 

(iv) The last term is 

+i?i(nM)Vap + 0(^2) = E{s)ip^p + ih^k{s)F{s, hu)ri{t)dMs, u) 
+hr](tfk{s)wlJo + 0{te), 
since R^{h,u) and Ri{h,u) are 0(^^) on the support of F. 

(v) The terms left combine themselves with the kinetic part and the time derivative of 
V'ap to give (|3.19j). 

Remark 3.3. The effective motion on M is given by the potential 

EJs) = (^ + V2) ^. X ^ ^ ^ X ^ ^ .3 23) 

a 

and is equal, also in this case, to the homogenized classical motion. 

The hypotheses that M is a nondegenerate critical curve and W is spectrally smooth 
imply that the normal oscillation is harmonic, and so all transversal states are equivalent. 
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3.3. The magnetic trap 

Using theorem ()3.H1 we can analyze the dynamics of a nonrelativistic particle in a strong 
magnetic field (magnetic trap). 

We suppose that the field is "strongly axially symmetric", i. e., that the vector 
potential is given, in cylindrical coordinates, by 

A{r,z) = £^{r,z)e. (3.24) 

The Hamiltonian is 



Since divA = 0, in the susbspace with zero angular momentum in the z direction 
fl3.25|) becomes 

H' = -T^-dArd^) - + ;r^^(r, z)\ 

2m r 2m 2mc^ 

or, scaling the wave function by the isometry 

V : L2(R^ X R,rdrdz) L'^{R^ x R,drdz) 

V-ijj = (3.26) 

VH^V^ = -^dl - - —dl + ^i/(r, zf. (3.27) 

2m 8mr^ 2m ' 2mc^ ^ ' ^ ^ ^ 

If we put m = 1 and consider the case of large electric charge, c/e = ah, we get in 

the end 

which, except for the centrifugal term, is of the form p.ip . with W{r,z) = £/{r, z)'^ /2. 

Theorem ()3.H1 tells us that, if we consider an initial state localized away from the 
origin, the semiclassical motion is constrained along the curve ^ (r, z) = 0, with effective 
potential given by 

E{s) = ndl[^~{s,u)r/2}^ll = ^\d,^~is,0)\ = ^\B{s,0)\, (3.29) 
where B is the magnetic field strength. 



e 



2 



4. Constraints by scaling of coupling constant: general case 

When the submanifold M has dimension (and codimension) greater than one, the theory 
developed in foregoing sections has to be generalized essentially in two aspects. 

First, if dimM > 1, the metric Gm, induced by the Euclidean metric of M""'"™' on 
M, may not be trivial, so both the classical motion of the centre of the squeezed state 
and the evolution of the dispersion matrices A and B have to be modified to take this 
into account. Thinking about the results we got above, it is not difficult to derive the 
new classical equations; we will simply obtain a motion on a Riemannian manifold with 
metric Gm{x) in the presence of a potential E{x) which is an eigenvalue of the reduced 
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Hamiltonian. In local coordinates this means (see, for instance, Abraham and Marsden 
1978) 

m = vit) (4.1) 

Vit) = -r{amv,v)-G~J{ait))VMa{t)), (4.2) 

where T{ri,riy = V^-jj]^rj'' {V'^-f^ are the Christoffel symbols associated to Gm) and V^, 
denotes the column vector whose coordinates are ft := d^ . 

The equations for the dispersion matrices are a bit more complicated, but, as we 
will see below, they can be derived, using the Hagedorn multiple scale method, from 
the term of order h of the formal expansion of the solution of the Schrodinger equation 
in powers of h^/"^. 

The second point is that, if codimM > 1, the Euclidean metric written in tubular 
coordinates is not diagonal anymore. 

In a formal expansion of the Hamiltonian if^ around the constraint, the off- 
diagonal terms give rise, as first noted by Maraner and Destri (Maraner and Destri 
1993; see also Froese and Herbst 2001, Mitchell 2001, Schuster and Jaffe 2003 and 
references therein) to an induced gauge field which minimally couples the longitudinal 
and the transversal motion. This gauge field is linked to the normal connection for 
the embedding M C R""*"'" (see, for instance, Spivak 1979) and it certainly vanishes if 
codimM = 1. 

At first sight, it might seem that in this case we can no longer split the motion into 
a tangential and a normal part, even in the semiclassical limit. 

Actually this is not true, since, applied to a squeezed state, the gauge coupling is 
of order and, due to the antisymmetric character of the normal fundamental form, it 
maps an eigenstate of the reduced Hamiltonian into a state which is ortogonal to it. 

According to the proof of theorem ()3.1|1 . this means that, if we start from an initial 
state which is concentrated along a classical trajectory, and we study its evolution when 
h goes to zero, the gauge term contributes only to the remainder and not to the leading 
term of the expansion in powers of /i^/^, which is again given by a wave packet in the 
longitudinal variables times an eigenstate of the normal Hamiltonian. 

In principle, higher order corrections can be calculated following the procedure 
developed by Hagedorn (Hagedorn 1994), even though in the general case the formulae 
can be cumbersome. 

In the following we will give some details of the calculations that justify these 
claims, even though, given the previous warnings, they are analogous to those of the 
two-dimensional case. 
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4.I. The Hamiltonian in tubular coordinates 

By the tubular neighbourhood theorem (Lang 1995), given a local chart : E C M 
M" for the submanifold M, and a 6 small enough, there exists a diffeomorphism between 

^(5) := {q e M"+'" : d{q, E) < 5}, (4.3) 

d{q,E) :=inf{|g-e| : e e E} 
and the open subset of the normal bundle of M given by 

TEj- := {(e, n) : e e n G T,E^ , \n\ < 5}. (4.4) 
The diffeomorphism can be chosen to be 

/(e,n) = e + ra, 

where we have identified every fibre T^E-^ with a subspace of M""*""^. 

This means that, given a (local) basis for the normal bundle {?^A;(e)}^i, we can 
write every point in f?(5) as 

q = ax) + yunk{ax)) x G M" y G (4.5) 

(summation over repeated indices is understood). 

Starting from the above expression, we can calculate the coordinate form of the 
basis for the tangent space in a point of ^{5) simply differentiating with respect to a 
coordinate Xi or yk, and then calculate the scalar product of two basis elements to get 
the local form of the metric. 

The result is 

r^f ^ f I n \ f GMiI~Sf \ f I nV 

^^^'^)=(o / j( /j(o / j ' ^^-'^ 

where 

iV,,,(x, y) = yuPt^x) /5f = ■ 9,n,, (4.7) 
Si,j{x, y) = yk{G]^l)iiafj{x) afj.(x) = Uk ■ ditj, (4.8) 

and ti denotes the basis for the tangent space (the index k and h always refer to the 
normal coordinates, while the other indices refer to the tangential coordinates). 

Pf^ and afi are called, respectively, normal fundamental form and second 
fundamental form of the submanifold M. Together with the metric Gm, they 
characterize completely the embedding of M into R'^"'""', up to a Euclidean motion 
(Spivak 1979). It is important to stress that = —Pi^, so, when codimM = 1, /? is 
identically zero. 

Using ()4.(j|l . we can write the Hamiltonian in tubular coordinates, but, as we did 
in the two-dimensional case, we have to modify the volume form given by g{x,yy^'^ := 
[det G{x, y)]^^"^, in order to get wave functions which have the right normalization when 
integrated over the submanifold M. 

After this, we have to dilate the normal coordinates by h, in order to separate the 
reduced Hamiltonian from the longitudinal part. 
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This can be achieved by the unitary operator 



{V^){x,y) = h-"^/'^{x,y/h) (4.9) 

\9{x,y)J 

V : L\cff{6/h),gM{xy^^dxdy) L'^{S{5), gi.x.yf'^dxdy) 



where guix) := detGuix). 
The result in the end is 

Hbo = V^HkV 



- ^p,(x, yr/%i" ( VI - V^N^ix, y), ^ ^ ) ^If pf 



-1/4 
Pa 



[/-;i5(x,i/)]-2G^/(a;) \ / V.-N{x,y)Vy 
I J y h~^Vy 

+V{x + hy) + {ah)-^W{x + hy), (4.10) 



where 



p.(.,y) = «<^ . (4.11) 

When we further expand the equation ()4.10|) . the terms containing pnix, y) give rise 
to additive corrections which depend only on the second derivatives (or the square of 
the first derivatives) of In p^j. They are of order at least h"^. This can be understood if we 
note that S{x,y) is linear in y, the second derivatives with respect to x are multiplied 
by ff' and 

lnp;,(x, y) = In ^^^{Gm[I -hS{x,y)f} ^ ^ ^^^^^^^ _ ^^^^^ ^ ^ ^^^^^^ _ 

det Gm 

= -2hTT{S) -h^TiiS^) + 0{h^) . (4.12) 

Therefore, in the following, we will put p/j = 1 without other comments. 

Expanding the potentials V and W, we obtain the reduced Hamiltonian 

^(2)(x) = -^-Ay + ^y^H{x)y + V{x), (4.13) 

where H{x) is the matrix of the Hessian operator in the basis {nk{C{x))}. The hypothesis 
that W has a smooth spectral decomposition implies that we can choose the Uk to be 
eigenvectors of H, so we can write y^H{x)y = J2xkx ^a(^)i/a,A;a- 

We will see in next subsection that, as before, the higher order terms in the Taylor 
expansion must be included in the remainder. 



4-2. The approximate evolution 

To construct approximate solutions to the Schrordinger equation 

ihd.ij = Hbo^ (4.14) 

we use the same procedure outhned in previous sections, which is based on the multiple 
scale expansion developed by Hagedorn (Hagedorn 1994). The operator ()4.10p is not 
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of the standard form studied in the hterature, so we briefly explain the modifications 
needed to cope with this case. 

When all the terms have been spelled out, ()4.1()j) has the form of an elliptic 
differential operator in x and y, with coefficients which depend on x and y as well 
as h, plus the reduced Hamiltonian, plus a remainder of order h, which comes from the 
Taylor expansion of V{x + hy) and W{x + %) up to first and third order, respectively. 

According to the Hagedorn method, to split the adiabatic and the semiclassical 
effects, we have to introduce a fictitious new variable 
X — ait) 

^ ■■= (4-15) 

which measures the "deviation" of the quantum evolution from the classical one, and 
consider ^ as an independent variable in the formal manipulations. 

Associated to ^, there is an auxiliary wave function, ip{x^ y,^;t), which satisfies the 
equation obtained substituting 

X — a{t) 



~( X — ait) \ 



into ()4.14j) . and adding to the right-hand side the term E(a{t) + h}^"^^) — E{x), which 
formally equals zero when = [x — a{t)]/h}^'^, where E{x) is a fixed eigenvalue of l)^'^\x), 
with multiplicity 1. 

When we perform this substitution, we replace the x dependence in the coefficients 
of the differential terms with a dependence on a{t) + h}^'^^. 

This is justified because when we apply a function of x, f{x), to a squeezed state 
(Pk{A, B, a, r], h, x), we can develop f{x) in Taylor series, up to order Z, around the centre 
of the packet, getting a remainder which, in norm, is of order (Hagedorn 1994 

and references therein). 

At this point, we make the Ansatz that 



i^{x,y,^]t) = exp [iS{t)/h) exp 



ivitfGMiam 



F{x, hy) 



y<9M{a{t))-'/\tPo + h'/^, + h^2 + .. .), (4.16) 
where a(t) and ri{t) satisfy equation ()4.1|) . S(t) is the associated action 

S{t) = ds^r]{sfGM{a{s))7]{s) - E{a{s)) , (4.17) 

and F is a smooth function which has support in x near the classical trajectory, and in 
hy near 0. 

Substituting this Ansatz in the equation for ip, and keeping terms up to order h, 
we can determine ipQ and , which, as shown in (Hagedorn 1994) are what is needed 
to solve the Schrodinger equation to lowest order in h^^"^. The calculations are lengthy 
and not very interesting, so we give simply the result. 

The approximate solution, up to order h^^'^, of ()4.14j] is 



^ap(a;, y; t) = exp (iS{t)/h^ exp 



gM[a 



{t))-^'^h-'''^ 



Constrained quantum systems 



17 



x^k{A{t), B{t), 1, 0, 0, y) + r{^) p^^V.$ 

+i'n{tfN{a{t),y)Vy^ + r^{tfGM{a{t))S{a{t), y)r^{t)^ 



|p|=3 



(4.18) 



where ^ is given by ()4.15p . and S are defined in ()4.7p and ()4.8p . and ^(x, is a real 
eigenstate of ^^'^\x), with eigenvalue E{x) of multiplicity 1. 

As before, f(x) is the bounded inverse of the restriction of [t)*^^^(x) — E{x)] to the 
orthogonal complement of in L'^iW^^dy). 

Remark 4.1. The evolution of the dispersion matrices A[t) and B{t) can be read from 
the terms of order h in the expansion, and contains explicitly the metric Gm'- 

d,A{t)u = r]k{t)[GMdpi}]u{a{t))A{t),i + i[Gl}{a{t))B{t)]a (4.19) 
d,B{t)a = '-r]{tf[GM{diG];l)GMMt)Mt)A{t),i + d^^E{a{t))A{t),i 

-Vk{t)[GMd,Gl^]kMt))B{t)ji (4.20) 
Remark 4.2. The term coming from the gauge coupling 

ir]{tfN{a{t),y)Vy<l> 

can be written, using creation and destruction operators for the normal oscillations, as 



{\,kx){v,K) 



d 



{\,kx){v,K) 



0J„ 



1/2 



- ^U.^UJ"^ • (4-21) 

Since /? is antisymmetric in (A, k\), (z/, h,^), the above expression is orthogonal to <1>, 
as we claimed in the introduction to this section. 



5. Takens chaos in quantum mechanics 



When the constraining potential is not spectrally smooth, that is, roughly speaking, 
when the eigenvalues or the eigenfunctions of its Hessian are not smooth, the classical 
motion on the submanifold M shows peculiar features. 

In this section we consider the quantum analogue of an example given by Takens 
(Takens 1980, see also Bornemann 1998) where W fails to constrain spectrally smooth. 

The Hamiltonian we study is 



+ 



pI, + P^ 



1 



(5.1) 



where q = {x,y) G M^, < ■, ■ > is the standard scalar product in M? and R{x) is the 
Rellich matrix (Kato 1995 and references therein) 



R{x) 



I 



Xi X2 
X2 -Xi 



(5.2) 
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The eigenvalues of R{x) are 

a;±(x)2 = ^(l±|x|), (5.3) 
with corresponding eigenvectors 

'"^'^ " 1^ sin(0/2) ) - 1^ cos(0/2) ; ' 

where = tan~^(x2/a;i), and the branch of the inverse tangent is chosen so that 
-7r/2 < < 37r/2. 

The eigenvectors are discontinuous along the semiaxis {x : xi = 0,X2 < 0}, or 
better, they exchange place upon crossing the cut. 

5.1. A brief review of the classical case 

To get a confining potential which is bounded from below we restrict the configuration 
space to 

S:={(x,i/):|x|<l/2}. (5.5) 
With this choice, the Hamiltonian (j5.1|) constrains the system to the submanifold 

M:={(x,y)GS:y = 0}. (5.6) 
An (almost) complete description of the limit motions when £ ^ is given by 
Theorem 5.1 (Takens 1980, theorem 3). Let 

W{q) = i < R{x)y,y>, 

then the solutions of the equations of motion 

1 

\/W(nJf.]] 

(5.7) 



Ut) = --VWiqeit)) 



qe{0) = g,(0) V. 

which satisfy 

Qv, ^ 0, (5.8) 

where Q : ^ M.'^ is the orthogonal projector Q{x,y) = x, converge uniformly to the 
unique solution of 

i(t) = -Vf/hom(a:(t),t) 
x(0) = x(0) = Qv„ 

where 

f/hom(a;,t) := t9+(t)cu+(a;) + 79_(t)cj_(x). (5.10) 

The functions 'd± are constant for t 7^ and can have any discontinuity in t = 0, 
provided that 19+ + remains constant. 

Remark 5.1. If Qv^ = 0, the limiting behaviour is not known. 
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5.2. A quantum analogue 



Hn = - — (A, + A^) + < g{\x\)R{x)y, y >, (5.11) 



In the quantum case, we consider the Hamiltonian 

f (A. + A.) + 

where g e C^(]R), g{z) = 1 when \z\ < 1/2, g{z) = when \z\ > 3/5. 

We use the same squeezing factor a for both transversal directions so that the 
eigenvalues of R{x) keep their simple form ()5.3|1 . 

The quadratic form < g{\x\)R{x)y,y > is non-negative, so Hf^ is essentially self- 
adjoint on Co°°(M4). 

Scaling y as we did in the above sections, we get 

Hbo = -^^x + i)ix) 

1 1 (5-12) 

^(x) = -^^y + 9{\x\)R{x)y,y> . 

Let us suppose from now on that \x\ < 1/2, so that = 1 (note that, in 

theorem (j2.ip . it is required that i){x) has an eigenvalue on an open set only, so this 
restriction is immaterial). 

To calculate the spectrum of f)(x) we exploit the fact that, for every x, R{x) is a 
real symmetric matrix, and can be diagonalized by an orthogonal transformation whose 
form can be derived from ()5.4|) . and is given by 

/oosW2) (5.13) 

^ ' y sin(0/2) cos(0/2) J ^ ^ 

It shows the same discontinuity of v±, but however is defined for all x. 
The corresponding unitary operator 



-1„ 



Z(x):L^(M^)-L 



[Z{x)^]{y) = ^iP{Z{x)-'y) 



(5.14) 



turns i){x) into the Hamiltonian of two uncoupled harmonic oscillators 



Z(i)tS(i)Z(i) = --A„ + —^^xfyl + ^_(x)V2- (5.15) 



The eigenvalues of [)(x) are then 

En+,n.{x) = Eofi{x) + —UJ+{x) + —UJ_{x) 

a a 

Eo,o{x) = = ^[(1 + \x\f^ + (1 - \x\)n 

5.2.1. The ground state The eigenfunction corresponding to £"0,0(2;) is 

$o,o(x,i/) = [Z{x)^>ofl]{x,y) = ^!Q,Q{x,Z{x)-^y), 
where \E'o,o is the eigenfunction of ()5.15j) with the same eigenvalue. 



(5.16) 
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The result, with a suitable choice of normalization constants, is 

c<j+(x)c<j„(x)"' 



exp(-^<i?(x)i/\y>). (5.17) 



The equations ()5.16p and ()5.17p tell us that both the energy and the wave function 
of the ground state of are C°° functions of x for \x\ < 1/2. Therefore, theorem 
()2.1|1 can be used also in this case, and gives us a constrained motion in the cylinder 
{{x,y) : |a;| < 1/2, y = 0}, with effective potential Eq^^x). 

The classical trajectory we obtained is the only one which is associated, in the 
funnel described by fl5.10|) . to a smooth homogenized potential. The semiclassical limit 
thus singles out a specific motion, which is linked to the initial normal oscillation. 

5.2.2. The excited states If we consider the excited states of f}(x), we observe crossings 
between different eigenvalues in x = 0. Unlike what happens in the classical case, 
however, an incoming semiclassical wave packet splits into two components only, giving 
rise to a bifurcation of the motion, and not to a funnel. 
For the first two excited states, for example, we have 

Eo,i{x) = Eo^oix) + uj.{x)/a (5.18) 
Ei^oix) = Eo^oix) + LJ+{x)/a, (5.19) 

(when |x| < 1/2 we have uj^{x) < 2uj-{x), so the other eigenvalues remain separated 
from these). 

The corresponding eigenfunctions are 

<l>o,i(x, y) = a-'/^%,o{x, y)[2uj4x)Y/^[- sin(0/2)yi + cos(0/2)2/2] (5.20) 
<l>i,o(x,y) = a-i/2^o,o(x,i/)[2cu+(x)]i/2[cos(0/2)|/i + sin(0/2)y2]. (5.21) 

Clearly, the two eigenvalues coincide when x = 0, and are not differentiable in such 
point, while the eigenfunctions are not even continuous. 

Carrying out a rotation between $o,i and ^i^, we can construct a smooth basis in 
the two-dimensional subspace generated by them. 

It is easily seen that 



<^A{x,y) \^ sin(0/2) - cos(0/2) \ / $o,i(x,y) 
Mx,y) J ■ \ cos(</./2) sin(0/2) J y <^i,o{x,y) 

are smooth in the origin, since [2c<j-|-(x)]^/^ = (1 ± Ixl)"*"/^ = 1 ± ^\x\ + 0(|x 



(5.22) 



2^ 



SO 



^A{x,y) = a ^/2<l>o,o(a;,y)| -yi- ^ViXi - ^y2X2 + Odx^j 
^B{x,y) = a"^/2<l>o,o(a;,i/)|i/2 + ^yiX2 - ^y2Xi + 0(|xn|. 



Note that 

< <!>B{x,y)Mx)<!>A{x,y) >l2(r2)= a"^ sin(</)/2) cos((/)/2)[cj_(x) - ^+(a;)] (5.23) 



a-' 



-^X2 + 0{\xf) ^0 Wx^O. (5.24) 
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Therefore, in Hagedorn's classification (Hagedorn 1994), this is a crossing of type 
/. The theory developed by him allows to elaborate on the qualitative features of the 
propagation we mentioned above. 

If the system is initially in a semiclassical state associated to the level Eq^i and 
passes through the region of crossing, a; = 0, with a non-zero velocity (this assumption 
of generic crossing was already present in Takens' theorem, ()5.8|) ) the final state is a 
superposition of two components, one evolving with the potential Eq^i and the other 
with the potential -Ei,o- More precisely we have 

Theorem 5.2 (Hagedorn 1994, theorem 6.3). There is an approximate solution 
'^{h,x,y,t) to the Schrddinger equation generated by the Hamiltonian i5.1^) that 
satisfies 

v[/(n,x,i/,t) = <l>o,i(x,y)exp(z5(°'i)^-(t)/n)^fc(A(°'i)^-(t),i?(°'^ 

+0{h}/^) (5.25) 
for t G [—T,Ti], for any Ti > 0. For t G [Ti,T], this solution satisfies 
<il{h,x,y,t) = <l>o,i(x,y)exp {tS'-''''>-'+ it)/h) 

X d^^''^^rn{A^'''^-^Ht), B^^'^y^^t), K a(°'^) (tW'^^ {t),x) 

m 

+<l>i,o(x,y)exp(5(i'°)'+(t)/;i) 

X E rf(;^'°V™(^^^'°)^+W,5(^'°)^+(t),^«^^'°Ht),^^^'°Ht),^) 

|m|<|fc| 

+0(^/2), (5.26) 
for some a > 0. 
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Appendix A. Nondegenerate critical submanifolds 

Let W : M"+™ ^ R be a non-negative function, and let M = {g G M"+'^ : W{q) = 0} 
be a smoothly embedded n— dimensional submanifold such that 

• M = {qe : DW{q) = 0}; 

• the Hessian H of W, defined as a field of linear operators H : M =Sf (/?""'"''") by 

< H{q)u, V >= D^W{q){u, v) u,v e g G M (A.l) 

(<■,■> is the standard scalar product in ]R"+™') is uniformly positive definite when 
restricted to TgM-^. 

Then, M will be called a nondegenerate critical submanifold of R"+"* and W will 
be called constraining to M. 
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Appendix B. Spectrally smooth constraining potentials 

Let W he a. potential constraining to a nondegenerate critical submanifold M. If the 
Hessian H oiW has a smooth spectral decomposition on M, 

r 

H{q)^J2'^,iqfP,iq), q e M, (B.l) 

k=l 

W will be called a spectrally smooth constraining potential. Here, uj^. and represent 
the (non-zero) eigenvalues and eigenprojections of the Hessian. 
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